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We explore the momentum-sector-selective metal insulator transitions recently found in the eight 
- site dynamical cluster approximation to the two-dimensional Hubbard model. The phase diagram 
in the space of interaction and second-neighbor hopping is established. The initial transitions from 
Fermi-liquid like to sector-selective phases are found to be of second order, caused by the continuous 
opening of an energy gap whereas the other transitions are found to be of first order. In the sector- 
selective phase the Fermi surface regions which are not gapped are found to have a non- Fermi-liquid 
self-energy. We demonstrate that the phenomenon is not caused by the Van Hove divergence in 
the density of states. The sector-selective and insulating phases are characterized by a cluster spin 
correlation function that is strongly peaked at the commensurate antiferromagnetic wave vector 
(tt, tt) but the model has no nematic instability. Comparison to dynamical mean-field studies on 
smaller clusters is made. 

PACS numbers: 71.10.Fd, 74.72.-h, 71.27.+a, 71.30.-|-h 



I. INTRODUCTION 

The "pseudogap" , a suppression of the electronic spec- 
tral function occurring for momentum states along the 
Brillouin-zone face but not for states along the zone di- 
agonal, is a basic and still ill-understood feature of hole- 
doped high-temperature cuprate superconducting ma- 
terials. (In electron-doped cuprates a different effect, 
confusingly also sometimes termed as pseudogap, is at- 
tributed to the presence of or proximity to long-ranged 
two-sublattice antiferromagnetic order.) The pseudogap, 
which occurs in the absence of any obvious long-ranged 
order, is a dramatic example of the more general phe- 
nomenon of "momentum-space differentiation," an in- 
crease in the variation in physical quantities around the 
Fermi surface as the insulating phase is approached. Its 
origin and consequences remain hotly debated topics. 

The cluster dynamical mean-field method^ pro- 
vides theoretical access to the momentum depen- 
dence of the electron self-energy and spectral func- 
tion. Application of this method to the two- 
dimensional Hubbard model has revealed strong indica- 
tions of momentum-space differentiation and pseudogap 
formationS^-i^^i^^^^^i^^^ as well as avoided!^ 
or actualii quantum criticality. In a previous paper— 
we demonstrated that when applied to an eight-site clus- 
ter the method yields a multistage, momentum-sector- 
selective metal-insulator transition in which different re- 
gions of the Fermi surface are successively gapped as car- 
rier concentration or interaction strength are varied. The 
successive gapping bears an intriguing similarity to the 
behavior of high- Tc cuprates in the pseudogap regime. 
In this paper we present a detailed analysis of this im- 
portant phenomenon. We determine the nature of the 
transitions, showing, in particular, that the initial tran- 




FIG. 1: Left panel: Brillouin-zone partitioning associated 
with the eight-site cluster including definition of the four in- 
equivalent momentum sectors A, B, C , and D. The non- 
interacting Fermi surface for t' = — 0.15f and density n = 1 is 
indicated by the gray line. Right panel: sketch of the para- 
magnetic state DCA phase diagram of the Hubbard model, 
calculated for the cluster shown in the left panel at half fill- 
ing, as a function of interaction strength U and next-nearest- 
neighbor hopping t' . A Fermi- liquid metal phase (left, red 
online), a sector-selective intermediate phase (middle, green 
online) and a fully gapped insulating phase (right, blue on- 
line) are shown. 



sition from a Fermi-liquid like phase to a sector-selective 
phase corresponds to the continuous opening of a gap (or 
pseudogap, for a discussion see Sec. lIXp in one momen- 
tum sector. The ungapped momentum sector is shown to 
exhibit non-Fermi-liquid features including a self-energy 
which vanishes more slowly than linearly as the frequency 
tends to zero. We determine how the behavior changes 
with carrier concentration, interaction strength, and vari- 
ations in band structure, clarify the effect of breaking 
particle-hole symmetry and show that the Van Hove sin- 
gularity does not cause the effects. We also present data 
indicating that the phenomenon is linked to a magnetic 
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instability and that the model is not unstable to nematic 
order. 

Most previous cluster dynamical mean-field studies 
have used interpolation methods to infer the behavior of 
lattice quantities from the limited information provided 
by the dynamical mean-field methods and analytical con- 
tinuation methods to infer the frequency dependence. We 
believe that it is essential to understand the behavior of 
the model directly; we therefore avoid "periodization" 
and continuation methods in this paper, basing our con- 
clusions on the analysis of directly measured quantities. 

As first noted by Biermann et alM> the phenomenon 
of momentum-space differentiation is related to the or- 
bitally selective Mott transition^!! which is found in the 
single-site dynamical mean-field theory of models with 
orbital degeneracy and has been extensively studied in 
the context of transition-metal oxides and actinides with 
degenerate d or /-levels as models for local-moment 
format ioni^i'^'^^ and as possible explanations of certain 
aspects of heavy fermion behavior.— The connection was 
explored by Ferrero et alr^ in a study of a two-site dy- 
namical mean-field approximation with a cleverly cho- 
sen cluster geometry and was also adopted by Liebsch 
et al^ In this paper we attempt to relate our results to 
the general understanding of the orbital-selective Mott 
phenomenon and make some comparison to the previous 
work of Ferrero et al. and Liebsch et al. 

An important question which we are not fully able to 
resolve is the relation between the sector-selective transi- 
tion we find and the onset of long-ranged order. We find 
an association between the gap and a pinning of density 
at a commensurate value which leads us to identify the 
transitions as Mott transitions. Our calculations are car- 
ried out in the paramagnetic phase of the model, where 
by construction no long-ranged order is possible. It is 
however quite conceivable that the onset of some sort of 
long-ranged order underlies the behavior we study. In 
four-site clusters this is believed to be the case: the gap- 
ping effect was related to the onset of short-ranged order 
of the plaquette singlet type in Ref. [13 while in Ref. a 
susceptibility analysis was used to argue that a transition 
to a state with long-ranged columnar dimer order domi- 
nated the physics of the four-site cluster. We present evi- 
dence here that antiferromagnetic correlations with range 
larger than the cluster size are important and we note in 
passing that nematic order appears not to be favored in 
the model we study. 

The rest of this paper is organized as follows. In Sec.HIl 
we define the model and crucial parameters and out- 
line the theoretical methods. In Sec. IIIII we explain how 
the phase boundaries are determined and in Sec. IIVI we 
present the calculated phase diagram. In Sees. |Vl IVIl 
and IVIII we discuss the physical content of our results. 
In Sec. IVIIII we place our results in the context of other 
work on orbital-selective transitions. Section HXl is a con- 
clusion, and is organized so that readers may turn di- 
rectly to this section, which summarizes the results with 
pointers back to appropriate portions of the main text. 



II. MODEL AND METHOD 

We study the two-dimensional one-band Hubbard 
model with Hamiltonian 

H ^^{fp- fJ-) cl^^Cp^a -f C/^ ni,|nj,4, (1) 

where 

Cp — —2t [cos{px) + cos(py)] — At' cos{px) coa{py). (2) 

The low {uj < 4eV)-energy physics of the high- 
temperature superconductors is believed to be described 
by a model of this sort with interaction U ^ 9t and next- 
nearest- neighbor hopping —0.3t ^ t' < —0.15t. A non- 
vanishing t' breaks particle-hole symmetry and we shall 
see that this has an important effect on the results. The 
carrier concentration n is controlled by the chemical po- 
tential fi and we shall be interested in dopings x ~ I — n 
between ±0.3. 

We use the "Dynamical Cluster Approximation" 
(DCA) formulation of cluster dynamical mean-field the- 
ory (Refs.[26land[ll) in which the Brillouin-zone is divided 
into N "patches" defined by the basis functions 4>a{p) 
which are 1 for p in patch a and zero otherwise and the 
electron self-energy is approximated as 

N 

LU) (j)a{p)^a{uj)- (3) 

Q = l 

The patches are required to have equal area and to 
fully tile the Brillouin-zone. Different DCA approxi- 
mations are defined by the number of patches and the 
choice of patch shapes and positions. We study here 
the N = 8-site DCA approximation defined by the tiling 
shown in Fig. [T] We label the patches by the 8 momenta 
^^^1=1.. .8 a-t the patch centers and distinguish inequiva- 
lent momentum sectors by letters A, B, C, and D as 
shown. An advantage of this tiling is that for the car- 
rier concentrations of interest the Fermi surface of the 
non-interacting model passes through two inequivalent 
sectors (labeled B and C). The cluster therefore gives 
direct (if coarse grained) access to momentum variation 
around the Fermi surface. In previous studies of two- 
sitoi^ and four-site^'i'^'^'^'^'^ i i ^ ^ i clusters, variations 
around the Fermi surface had to be inferred from data 
involving also the momentum sectors (0,0) and (tTjTt), 
which are far from the Fermi surface so that variation 
away from the Fermi surface was mixed in with variation 
around the Fermi surface. 

The frequency-dependent functions £0(0;), a = 
1 , . . . , 8 are obtained from the solution of an eight- 
site quantum impurity model^i^l defined by an action 
S ~ J drdr' C{t, t') which is a function of a continu- 
ous (imaginary)-time variable and the momentum sector 
labels Ki=i...s corresponding to the centers of the tiles 
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shown in Fig. [TJ Specifically, 

^ = J2^Kir^r')d^i,Jr)dK,Ar') 



(4) 



K 



Ki...Ki 

with ^ denoting a sum over all K such that Ki + K2 + 
+ equals a reciprocal-lattice vector. The mean- 
field functions are obtained from 



{dk)G{k) 



K 



(5) 



with Jj^{dk) standing for an integral over the tile cen- 
tered on momentum K, normalized so that ^^{dk)\ ~ 1 
and G{k) is the Green's function of the lattice problem 
computed with S defined by Eq. jS]). Note that within 
the DCA approximation the different momentum sectors 
are not mixed by the self-consistency condition. 

Following the suggestion of Ref . [l^, Ferrero et ali^ and 
Liebsch et al^ observed that the DCA equations have 
the same form as those used in the single-site dynami- 
cal mean-field theory of a multiorbital system such as a 
transition-metal oxide or the / levels in a heavy fermion 
material. In the latter case one identifies the labels K 
with the different local orbital states, the Jj^ becomes 
the integral over the whole zone of an appropriate pro- 
jection of a multiband lattice Green's function, and the 
self-consistency equation is in general not diagonal in the 
indices K. 

We solve the impurity model with the numerically ex- 
act continuous-time auxiliary field technique^^ with de- 
layed updatesi^ Calculations are rendered difficult by 
a fermion sign problem which causes the computation 
to scale exponentially instead of polynomially in inverse 
temperature. The sign problem vanishes in the particle- 
hole symmetric case t' = 0, n = 1 but becomes more 
severe as t' and |n — 1| are increased. The severity of the 
problem may be characterized by the average sign (s) as- 
sociated with the simulation of a fermionic problem (see, 
e.g., Ref. [30! ) because errors scale roughly as l/(s). In 
sign-problem-free situations (s) = 1 at all temperatures 
but in the typical cases studied here (s) < 1 and decreases 
exponentially as temperature T —^ 0. The severity of the 
sign problem increases with U and <' and is apparently 
more severe for moderate than for small or large doping. 
For example, at interaction strength U = It, t' /t = —0.3 
and inverse temperature jSt = 25 wc find (s) ~ 0.15 for 
fi/t = 1 (around 5% hole doping, near the boundary of 
the sector-selective transition) but (s) ~ 0.85 at U/t = 8, 
fi/t = 2 far from the boundary. The sign problem is also 
strongly dependent on the cluster geometry. 

The lowest temperature accessible with the computa- 
tional resources available to us is /3< « 40 at C/ < 8t and 
most of our results are obtained at f3t ~ 20. The con- 
ventional high-Tc band parametrizationSi implies that 
Pt = 20 corresponds to temperature T 200K. 
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FIG. 2: Total density per spin as a function of chemical 
potential, for U/t = 7, pt — 20, and second-neighbor hop- 
ping t'/t — —0.15 (solid line and circles, black online) and 
t' /t = —0.3 (dashed line and diamonds, red online) 



Throughout the paper we mainly present results as a 
function of chemical potential. For reference, we show in 
Fig. [2] the total density as a function of chemical poten- 
tial, ior U/t = 7, pt = 20, and two representative values 
of the next-neighbor hopping t' . 



III. ANALYSIS 

As noted above we base the conclusions of this pa- 
per on the analysis of directly measured quantities. 
The most important and easily accessible of these are 
the imaginary-time sector Green's functions Gk = 
(TtCk{t)cj^(0)^ = {Gx^ — ^k) ^, from which we also 
may extract the sector density as 



UK 



-GA-(r^O-). 



(6) 



Note that the total particle density is n = tik/S. 

It is useful to consider the value of the sector Green's 
function at the center of the imaginary-time interval, 
Gk{t = §)■ Use of Eq. ^ and the usual spectral repre- 
sentation formulas shows that 



0Gk ('f 



duj 



2t:T cosh[a;/(2T)] 
with sector spectral function 



(7) 



Ak{lu) 



{dk) — 

K 



ImE/f (w) 



fi-Ek- Rel]/f(w))2 



ImS|.(w)' 



(8) 

The integral in Eq. ([7]) is dominated by a; < T so that 
in the low-temperature limit the measurement provides 
information about the Fermi-level sector-averaged elec- 
tron spectral function. We may distinguish three cases. 

In a Fcrmi-liquid metal, lmT.(uj s- 0) is negligible and 
/i— Efc— ReSif (w = 0) vanishes for some k in the sector. In 
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FIG. 3: Upper panel: ^G(%) for sectors B (full symbols) 
and C (empty symbols), at U jt — 7 and t' /t = —0.15. The 
strong temperature dependence in the sector C curves arises 
from the Van Hove divergence in the density of states. The 
crossing points indicate the onset of gapping in the sectors. 
Lower panel: same, for t' /t — —0.3. 



this case f3GK{P/'2) becomes the non- interacting density 
of states Anon = J[({dk)S{ficS — Sfc) evaluated at the 
renormalized chemical potential 



11°^ - 



/i-ReSA'(w = 0). 



(9) 



One important caveat is that sector C contains the Van 
Hove singularity of the two-dimensional band structure, 
where the density of states diverges. This divergence 
leads to a strong temperature dependence in f3G{f3/2) so 
that the non-interacting density of states value may be 
reached only at very low temperatures, often below our 
measurement limit. 

In a non-Fermi-liquid metal ImS(ti;) vanishes less 
rapidly than lu and fj, = Sk — ReEi<-(w = 0) is satis- 
fied for some k in the sector and as T — *■ (3Gk{P/'2) 
takes a value which is nonvanishing but typically less 

than Anon(Moff)- 

Finally, we may have a gapped state if ImE(Li; ^ 0) ^ 
and fj, — Sk — ReEif (w = 0) docs not vanish for any 
k in the sector so that f3GK{P/'2) vanishes exponentially 



at low T. We shall see that in our calculations this ef- 
fect arises from the appearance of a low-frequency pole in 
the self-energy which pushes states away from the Fermi 
level, making it difficult to satisfy the quasiparticle equa- 
tion at low energies. 

The boundaries of the gapped regions may, in princi- 
ple, be determined from any of the density, the Green's 
function, or the self-energy. In practice significant uncer- 
tainties, compounded by our limited temperature range, 
exist in each method. We believe that the most reliable 
method is to plot /3G(/3/2) as a function of interaction 
strength or chemical potential for several temperatures 
as shown in Fig. [31 

For sector C on the negative fi (hole-doped) side a 
clear crossing point at fi « ~l.5t is evident, separat- 
ing a regime where (3G{(3/2) increases as T is decreased 
from a region where f3G{P/2) decreases as T is decreased. 
We interpret the crossing point as marking the chemi- 
cal potential at which a gap begins to open in sector C. 
We observe that for < — 1.5i, /3G(/3/2) is substantially 
less than the Fermi-liquid value but increases as T is de- 
creased. This behavior is consistent with the hypothesis 
of a Fermi liquid state with very large thermal correc- 
tions due to the divergence in the density of states at the 
Van Hove point. At the temperatures available to us the 
possibility of a marginal Fermi-liquid or non- Fermi-liquid 
state in this parameter regime as proposed in Ref. |T3 can 
neither be ruled out nor confirmed. 

In sector B the identification of the transition point in 
the doping case is more complicated (in the interaction- 
driven case an analysis as in sector C is straightforward) . 
As can be seen from Fig. |3] (see also Fig. [TT|) a sim- 
ple crossing point does not occur in the (3G{(3/2) graph. 
Rather, one begins to see a fan out from a temperature- 
independent set of curves. Furthermore, in the range 
— 1.3i < /i < — 0.9i (3G{(3/2) seems to evolve (see also 
Fig. [TT|) as T is decreased to a value which is temper- 
ature independent but less than the Fermi-liquid value 
(3Go{(3/2, ficfi)- This is evidence for a non-Fermi-liquid 
state. As the chemical potential is further increased, ev- 
idence of a temperature-dependent decrease in /3G'(/3/2) 
becomes apparent, and by fi = —0.5t a clear gap has 
opened. Precisely locating the point at which the physics 
changes from a gapless non-Fermi-liquid to a gapped 
state is thus challenging but it is clear that the onset 
of gapped behavior in sector B occurs at a substantially 
higher chemical potential than the onset of such behavior 
in sector C. 

One may also consider density vs chemical-potential 
curves such as those shown in Fig. Interpretation is 
complicated by two issues: first, if the gap is small then, 
at the temperatures accessible to us, a substantial tem- 
perature variation occurs. Second, while our data indi- 
cate that if sector K has a gap, then at T = uk = 1/2, 
we see that the density may approach the pinned value 
from above or below depending on the value of the chem- 
ical potential. For these reasons we do not use the nx 
data to identify phase boundaries. 
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Alternatively, the phase boundaries can be determined 
by the evolution of the self-energy or inverse self-energy 
with temperature. We chose not to use this method, as 
self-energies in the sector-selective region display a strong 
temperature dependence that can be analyzed more ac- 
curately by considering the crossing point of /3G(/3/2). 



IV. PHASE DIAGRAM 
A. Interaction-driven transition, half filling 

Figure [5] shows the phase diagram determined at half 
filling as a function of interaction strength U and second- 
neighbor hopping t' (note that at n = 1 the phase dia- 
gram depends only on \t'\). For t' not too large, a two- 
stage transition occurs, in which as the interaction is in- 
creased a gap first opens in sector C followed by a sec- 
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FIG. 4: Upper panel: Sector-specific density n(/i) for U/t = 7, 
f3t = 10,15,20,25, and t' /t = -0.15. Open symbols: Sec- 
tor C. Filled symbols: sector B. Insets: expanded view of 
the hole-doped (upper left) and electron-doped (lower right) 
metal insulator transition regions. Lower panel: same for 
t' /t = —0.3. The pinning of the density to 1/2 in the insulat- 
ing phase is clearly visible. Arrow: location of the first-order 
transition. 
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FIG. 5: Phase diagram calculated at half filling in the plane of 
interaction strength U and second-neighbor hopping t' from 
measurements of the crossing point of j3G{(i/2) at jSt — 15 
and 20 as described in Sec. IIIII Lower region (open squares, 
red online): "metallic" region with gapless behavior in sec- 
tors B and C. Upper region (stars, blue online), "insulating" 
region with all sectors gapped. Shaded area (filled squares, 
black online): sector-selective region with sector C gapped 
and sector B gapless. 



ond transition to a completely gapped state. However, as 
the second-neighbor hopping amplitude is increased, the 
range in U over which only one sector is gapped decreases 
and for \t'\ > 0.25t there is only one transition. 

Wc may relate the orbitally selective behavior to the 
density of states ^{e) = J {dk)TTd{e — Sk) shown for the 
two relevant sectors in Fig. [S] At t' = the densities of 
states of sectors B and C are symmetrical about e = 



— t'/t = -0.3 

— t'/t = -0.15 

— t'/t = 




Energy 



FIG. 6: Noninteracting density of states computed from 
the analytically known /3G'o(/9/2, /i) at inverse temperature 
P = 25t for t'/t = 0,-0.15, and -0.3. Bold lines: sector 
B. Dashed lines: sector C. Sectors A and D are not shown. 
Arrows indicate the chemical potential corresponding to half 
fiUing. 
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(the chemical potential corresponding to half filling) but 
the density of states in sector C is sharply peaked at e cor- 
responding to the Van Hove energy £{px — 0,py = n/2) 
(= at i' = 0) and is narrower than that of sector B. It 
is natural that the sector with the narrower band and the 
higher density of states should undergo the Mott transi- 
tion first As t' is changed from zero the posi- 
tion of the Van Hove point shifts in energy, although the 
positions of the band edges of sector C do not change 
[the bandwidth of sector B does change from A\/2t — At' 
to 4i for \t'\ > {y/2 — The chemical potential corre- 
sponding to half filling similarly shifts. The existing lit- 
erature on the orbitally selective Mott phenomenon does 
not provide guidance on the fate of the sector-selective 
transition in these circumstances. We note, however, 
that a metric for estimating the location of the Mott 
transition is to compare the interaction energy to the 
U = Q kinetic energy given by J^^^^^^isp - fJ.)f{ep - fi) 
with / the Fermi function. As t' is shifted from zero 
to — 0.3i the kinetic energies of the B and C sectors 
change from (-0.084,-0.018) to (-0.035,-0.023); thus 
as t' is changed from zero the kinetic energies of the bands 
become similar, explaining the disappearance of the or- 
bitally selective phase. 



B. Doping-driven transition, fixed interaction 

Sector-selective transitions occur also as functions of 
doping at fixed interaction. For interaction strengths ly- 
ing in the orbitally selective (shaded) region of Fig. [5] 
sector B remains gapless as the chemical potential is 
changed whereas there is a critical doping level (or chem- 
ical potential) beyond which the sector C gap closes. We 
do not consider this case further, focusing instead on U 
large enough that both sectors are gapped at half fill- 
ing. Figure [7] presents a phase diagram in the space of 
chemical potential and t' for U = 7t. For small \t'\, both 
electron and hole doping are seen to produce a two-step 
transition, in which the first stages of doping take place 
in sector B and only subsequently does sector C start 
to be doped. Interestingly, as —t' is increased a qualita- 
tive particle-hole asymmetry develops. On the electron- 
doped side the two transitions coalesce (within our res- 
olution) to one and become first orde r ^'^i^^ while on the 
hole-doped side the presence of two transitions appears 
to be a generic feature. 

V. CHARACTERIZATION OF TRANSITIONS 
A. Interaction-driven transition 
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FIG. 7: Phase diagram at interaction strength U = 7t as & 
function of chemical potential ji/t and second- neighbor hop- 
ping t' /t, from measurements of l3G{(i/2). Positive /i corre- 
sponds to electron doping and negative /i to hole doping. A 
particle-hole transformation converts t' —t' so it is suffi- 
cient to show only one sign of t' but both polarities of doping. 
Dashed lines (circles and squares) indicate chemical poten- 
tials at which the C-sector gap opens, heavy solid lines (dia- 
monds and triangles) indicate chemical potentials where the 
B-sector gap opens. Note that for —0.15^ the critical chemical 
potential ^off at which the sector-selective transition occurs 
corresponds to a hole doping of x ~ 0.109. Thick black line: 
location of the coalesced first-order transition. 



In this section we characterize the interaction-driven 
transitions found in the eight-site cluster. We begin by 
summarizing what is known about metal-insulator tran- 
sitions in other dynamical mean-field theory implemen- 
tations. In the single-site dynamical mean-field approx- 
imation the correlation-driven metal-insulator transition 
occurs at a large U and has a complicated structure in 
which a preformed, large-magnitude gap exists on both 
sides of the transition and metallic vs insulating behav- 
ior is determined by the presence or absence of midgap 
statesi^ On the other hand, in four-site cluster DMFT 
approximations the transition occurs at a relatively small 
U ^ 5t and is apparently always first order, characterized 
by the discontinuous opening of a gap as the correlation 
strength is increased above a critical value.— 

Figure [S] shows our results for the sector density of 
states for several t' values. We observe, in agreement 
with Fig. 2 of Ref. [H which presents data for only two 
t' but at the lower temperature /3i = 40, that for t' near 
zero the two transitions are well separated in U and that, 
while the C-scctor transition appears smooth, marked 
by a crossing point in f3G{f3/2), the i?-sector transition 
seems discontinuous, marked by an apparent jump in 
G{(3/2). As the magnitude of t' increases, the C tran- 
sition steepens, and beyond the point where the sector- 
selective phase disappears the single transition becomes 
first order as in the four-site case. 

We may understand the change in transition order 
from consideration of the sector densities. The breaking 
of particle-hole symmetry caused by a nonzero t' means 
that in the weakly correlated metallic phase at total den- 
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FIG. 8: l3G(P/2) at density n = 1 for sector B (filled symbols) 
and C (open symbols) as function of interaction strength U/t 
for several values of the next-nearest-neighbor hopping t' /t at 
(3t = 15. At t'/t = -0.3 f3t = 20 data are also shown (thin 
line, small symbols). 



sity n = 1 ric ^ 1/2. For t' small, nc is close to 1/2 and 
evolves smoothly to the pinned value as U approaches the 
sector-selective value. However for larger t' the difference 
of nc from 1/2 is too large, and the evolution with U is 
preempted by a first-order transition. 

We now show how this behavior is reflected in the self- 
energy. We expect the opening of a gap to be related to 
the appearance of a pole in the self-energy so that 
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ImSA'(w) = 7rA25(cj ~ UJ^) + Kegi^) (10) 

implying 

^KiioJn)^-. j^+Ao+AiiuJn + ---, (11) 

ILUn — iiJp 

where Aq^i and the ellipsis arise from the regular part 
Sj-eg of ImE. Reference |35| found that correlation gaps in 
Hubbard-like models were reliably estimated from the so- 
lution to the quasiparticle equation cj — £^-1-//— ReS(ix') = 
0. If we neglect the A terms we may solve the quasiparti- 
cle equation to obtain an estimate for the sector half gap 
in terms of the band-edge energy ~ 2\/2t — 2t' 
[sector B, \t'\ < {■s/2 - \)t] and 2 (sector C), 

^{Eii + + 4A2 + ^{E^ - ujpY + 4A2 - 2E' 

E,^ 2 

(12) 

In the particle- hole symmetric case = = 0, so 
Yiiiujn) is purely imaginary and would diverge as ujn — > 
while Eg = ^(i;^')2 +4A2 - E'^ . Figure [S] reveals es- 
sentially this behavior. In sector C at U = It drawing a 
straight line through the origin at the lowest Matsubara 
point yields A2 w 4/^2, implying Eg 1.25t, reason- 
ably consistent with the gap defined from the chemical 



FIG. 9: Imaginary part of the inverse of the Matsubara-axis 
self-energy for momentum sectors B and C, measured as a 
function of interaction strength U /t at density n = 1, t'/t — 
and inverse temperature pt = 40 (/3i = 15 for U/t — 8). 
Solid lines indicate fits to U/t = 6.75, 7, and 8 with slopes 
ImE^^(LJo)/t^o = 0.14,0.25, and 0.30 (in units of f"^), re- 
spectively. 



potential in Fig. [7] while in sector B a similar analysis 
gives (with less confidence) 0.7t, slightly smaller than the 
chemical-potential gap. 

As the interaction strength is decreased from U ^ It 
to U ~ 6.55^ the C- and i?-axis self-energies and there- 
fore gaps change only slowly. However, as the interaction 
is further decreased from 6.55i to 6.45i to 6.4t the evi- 
dence for a gap in sector B vanishes abruptly and the 
gap magnitude in sector C becomes smaller. As the in- 
teraction strength is yet further decreased the C-sector 
gap drops rapidly. The frequency range implied by the 
temperatures shown does not permit a reliable gap analy- 
sis, and, in particular, does not exclude the possibility of 
a pseudogap such as discussed in Ref. [l^ but the rough 
estimate obtained by drawing a straight line from the 
origin through the lowest Matsubara frequency indicates 
that the gap (or pseudogap) at C/ = 6.2t is about a factor 
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of four smaller than that at J7 = 6.65t. 



-1 



-2 





1 

•-• Sector B 

♦— » Sector A 
Sector D 
■-■ Sector C 
- - Van Hove 




t 


- 


1 







_l , U I 1 L 

-4 -2 

ll/t 



FIG. 10: Sector-dependent effective chemical potential /ieff — 
/i — ReE(aJ = 0), plotted as a function of chemical potential for 
U/t = 7,t'/t = -0.15, and fit = 20. The horizontal dashed 
line indicates the ^eff value corresponding to the Van Hove 
point uvu ~ —0.6 in the C-sector density of states (Fig. |6]) 
and the vertical dashed line indicates the /j, value at which 
fJ'cS,c = MvH- The "effective" location of the Van Hove sin- 
gularity in sector C is shifted to fi/t = —2.1. Values in the 
insulating region are excluded. ReE was obtained by linearly 
extrapolating the first two frequencies. 



B. Doping-driven transition 

In this section we characterize the doping-driven tran- 
sition at a correlation strength sufRcicntly strong that 
both sectors are insulating at n = 1. We begin with an 
analysis of the dependence of density on chemical po- 
tential, shown in Fig. [D The first noteworthy point is 
that on the hole-doping side the sector C transition oc- 
curs at the chemical potential corresponding to sector 
C being half filled. The upper left inset shows an ex- 
panded view of the C-sector occupancy for different tem- 
peratures, demonstrating that deviations from the half- 
filled value are due to thermal excitations in a particle- 
hole asymmetric situation. We have found similar results 
for all other parameters and for this reason identify the 
sector-C transition as an orbitally selective Mott transi- 
tion. 

The important role played by the particle density is 
further demonstrated by the lower right inset of the up- 
per panel of Fig. |4l which provides an expanded view of 
the electron-doped transition. We see that as T — *■ the 
C-sector occupancy apparently develops a discontinuous 
jump. We observe that in the non-interacting model at 
t' < 0, half filling corresponds to the C sector being more 
than half filled and the B sector being less than half filled, 
with the filling difference increasing as the magnitude of 
t' increases. Thus as we approach half filling from the 



electron-doped side, encountering a sector-selective tran- 
sition means shifting electrons from sector C to sector 
B. For t' near this is possible; for larger amplitude t' 
this becomes too energetically expensive. The two tran- 
sitions coalesce into one first-order transition. This can 
be seen in the lower panel of Fig.[5]for t' /t = —0.3 which 
shows that the density of the B sector decreases below 
the value of half filling when the chemical potential is 
raised (indicated in the figure by an arrow) while the 
C-sector density jumps vertically to around 0.59. 
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FIG. 11: Zoom of Fig. O estimator l3G{l3/2) for the many- 
body density of states calculated as a function of chemical 
potential fj, for U/t = 7 and t' /t — —0.15 for sector B (upper 
panel) and sector C (lower panel). Solid line (purple online): 
Fermi-liquid value j3G{/3/2) evaluated for noninteracting elec- 
trons at the chemical potential /iotr. 



An alternative view of this physics is provided by 
Fig. 1101 which displays the renormalized chemical po- 
tential fics defined from the zero-frequency limit of the 
self-energy via Eq. ^ for the different sectors and t' = 
— 0.15<. For low carrier concentrations /x < — 3t the self- 
energy is almost momentum independent and the differ- 
ence between //off and fi simply reflects the correlation- 
induced renormalization of the n{ij,) curve familiar from 
elementary analyses. As fi is increased through — 3t 
momentum-space differentiation begins to occur. Of par- 
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ticular interest iiere is tiiat the C-sector ^^s lies above 
the B-sector /Xeg, signaling the onset of a relative shift of 
carriers from C to B. The dashed hnes indicate the renor- 
malized chemical potential at which the C-sector Fermi 
surface passes through the Van Hove point. Figure [10] 
shows that the ^esifJ-) curves pass through /icff = /ivH 
with no noticeable change in behavior. Figure [TT] shows 
that the many-body density of states evolves smoothly as 
/i is increased through the point where ^cft = ^J■YH■ The 
solid curve in the lower panel of Fig. [TT] shows (3G{(3/2) 
computed for non-interacting electrons at ^ = ^off- The 
Van Hove peak is evident and occurs at a significantly 
lower than the crossing point ^c" — — 1.5t, at which 
sector C becomes gapped. This physics is also signaled 
by the rapid upturn in ^oS,c^ which occurs because for 
/i > nc is pinned so the remaining doping is forced 
into sector B. Finally we note that on the electron-doped 
side the relative ordering of the ^off is opposite, signaling 
a transfer of carriers from C to B. 

The upper panel of Fig. [TT] reveals an additional inter- 
esting feature: a range from /i w — 1.5i to w — 0.7t, 
where the estimator for the _B-sector density of states, 
/3Gs(/3/2) apparently reaches a stable low-T limit which 
is neither zero nor the Fermi-liquid value. Examination 
of the fT.(/i) plots in Fig. [T] confirms that the sector B 
transition point is « — O.St so that in the whole non- 
Fermi- liquid range — 1.5t < /i < — 0.5t sector B is not 
gapped. We interpret this as a non-Fermi-liquid phase 
for sector B, which is approximately coincident with the 
sector-selective phase. The temperature dependence vis- 
ible in the sector C plot suggests that for < — 1.5t both 
sectors have a Fermi-liquid ground state, although much 
lower temperatures would be required to definitively es- 
tablish this. Non-Fermi-liquid physics in orbitally se- 
lective Mott regimes of models of orbitally deg enerate 
transition-metal oxides was reported by Refs. [l3 and [2^. 
However, our finding of more Fermi-liquid like behavior 
in the orbitally selective regime of the interaction-driven 
case suggests that there is not a generic association be- 
tween sector selectivity and a non-Fermi-liquid behavior 
of the ungapped sector. 

We are aware of three possible physical interpretations 
for the apparent non- Fermi-liquid phase indicated by this 
analysis. It may be intrisically non-Fermi-liquid, indicat- 
ing that over a nonzero range of the Fermi surface the 
imaginary part of the electron self-energy vanishes less 
rapidly than cj as w — > 0. Alternatively, it might be 
that the gap which is observed in sector C is also present 
over some of the momentum range included in sector B 
so that the measured non-Fermi-liquid behavior arises as 
an average over some parts of the Fermi surface which 
are gapped and some which are gapless and Fermi-liquid 
like. Finally, we can of course not rule out the possibility 
that there is a Fermi-liquid scale which is simply far be- 
low our measurement temperatures. Distinguishing these 
possibilities requires studies of larger clusters and lower 
temperatures, which are at or beyond the limits of our 
present computational capabilities. 
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FIG. 12: Real (upper panel) and imaginary (lower panel) 
parts of sector C self-energy calculated for indicated chemical 
potentials corresponding to hole doping, at U/t = 7, t' /t = 
—0.15, and pt = 20 (thick lines, thick symbols) and pt = 25 
(thin hnes, thin symbols). 



We now turn to the behavior of the self-energy. Fig- 
ure [T^] shows the real and imaginary parts of the (7-sector 
self-energies for the doping-drivcn transition. As is in- 
creased above nc ~ —l.Qt both real and imaginary parts 
develop the upturn expected if a pole appears in S. The 
temperature dependence given for ImS confirms this in- 
terpretation. Figure [13] shows the corresponding plots for 
sector B. While much lower temperatures (not accessi- 
ble with presently available computer resources) would 
be required to establish the precise behavior, it is clear 
that for a range of ^ around —1 ImS^ does not seem 
to linearly extrapolate to zero while neither ImS^ nor 
ReEs has a pole behavior. For /i < —IM Fermi-liquid 
behavior is apparently restored while for ^ > — 0.6t a 
pole begins to develop. 
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FIG. 13: Real (upper panel) and imaginary (lower panel) 
parts of sector B self-energy calculated for indicated chemical 
potentials corresponding to hole doping, a.t U /t = 7 , t' /t = 
—0.15, and [3t = 20 (thick lines, thick symbols) and pt = 25 
(thin lines, thin symbols). 
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FIG. 14: {SzkSzk') for K = K' = {n,n) at U/t = l,l3t = 15, 
and t' /t = —0.15 for the eight-site (bold lines, filled symbols) 
and four-site (dashed lines, empty symbols) cluster geome- 
tries, as a function of doping. 



frequency component grows substantially (by a factor of 
4) and the first Matsubara frequency grows somewhat 
(factor of w 2) while the other frequencies are essentially 
unchanged. Thus we see that the sector-selective regime 
is associated with a large value of quasiclassical (r2„ = 0) 
spin correlations. This finding is consistent with the ar- 
guments of Kyung et al^ that the pseudogap is associated 
with backscattering arising from slow, reasonably long- 
ranged spin fiuctuations. By contrast, the correlation 
functions of the four-site cluster have a weaker frequency 
and doping dependence, as expected from the dominance 
of the singlet correlations on the plaquetteJ^ However, 
it is not clear whether a spin fluctuations approach can 
account for the pinning of the density to the commensu- 
rate value of 1/2 that we observe to be a hallmark of the 
sector-selective phenomenon in the eight-site cluster. 



VI. SPIN CORRELATIONS 

To obtain further insight into the physics of the sector- 
selective transition we present in Fig. [14] the cluster 
spin-correlation functions obtained from S^(t — t') = 
{Sf{T)Sj{T')), where i and j are cluster site indices. 
While the cluster correlation functions are not identi- 
cal to the spin correlators of the lattice model, they are 
expected to show the same physics. We have Fourier 
transformed the correlation functions into Matsubara fre- 
quency and cluster momentum-space. For all cluster mo- 
mentum sectors except (tt, tt) the spin-spin correlation 
function is small (four times smaller than S^^.n at fx = —2t 
and 20 times smaller at fi — —O.M) and only weakly fre- 
quency and doping dependent. We therefore present in 
Fig. [14] only the (tt, tt) correlators. For comparison we 
show also the same correlator for a four-site cluster. We 
see that as chemical potential is decreased through the 
sector-selective value nc ^ —1.6t the zero Matsubara 



VII. ABSENCE OF POMERANCHUK EFFECT 

A nematic instability^ corresponding to a breaking of 
rotation symmetry without a corresponding long-ranged 
order has been extensively discussed in the context of 
strongly correlated systemSi ^'''^^i'^^'^°i'^^'^^i'^^i"^^ For ex- 
ample, in the context of the metamagnetic material 
Sr3Ru207 a transition corresponding to an increase in 
density of electrons with momentum near (tt, 0) and cor- 
responding decrease in the density of electrons with mo- 
mentum near (0, tt) (or the reverse) has been discussedi^ 
The locking of the sector density to the commensurate 
value nc = 1/2 in the sector-selective region motivates 
the question whether the nematic fluctuations are partic- 
ularly large near the sector-selective phase. A tendency 
toward nematicity would correspond to an anticorrela- 
tion of density fluctuations between different symmetry 
equivalent sectors. Figure [TJ] shows the results of our 
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FIG. 15: Absence of Pomeranchuk effect: density-density cor- 
relators for U/t — 7, pt = 15, and t' /t — —0.3 for momenta 
K 7^ A''. All correlators are zero within error bars. 



measurements of the appropriate correlation function as 
a function of chemical potential. We see (within errors) 
no correlations, even near the sector-selective point. We 
therefore conclude that within the eight-site DCA ap- 
proximation the Hubbard model does not exhibit any 
tendency toward charge nematic behavior at the inter- 
mediate couplings we have studied. This finding is in 
disagreement with RG studies of the Hubbard model^i^ 
and simulations of the t-J modelj^i^i^ where a Pomer- 
anchuk instability is seen. Whether the discrepancies are 
a consequence of a difference in parameter regimes exam- 
ined or arise from differences in approximation methods 
remains an open question. 



VIII. RELATION TO ORBITAL SELECTIVE 
MOTT TRANSITION 

In the DCA implementation of cluster dynamical 
mean-field theory, the impurity-model Green's function 
is diagonalized in a basis labeled by momentum quan- 
tum numbers K which are related, via a self-consistency 
condition, to averages over particular regions of the 
Brillouin-zone of the lattice model. The model may be 
viewed more abstractly, regarding the labels K as or- 
bitals of a general multi-level impurity model and aver- 
ages over different regions of the Brillouin-zone as arising 
from different baths. In this way, as noted by Bicrmann 
et al.f^ the DCA equations may be mapped onto the gen- 
eral problem of a multilevel impurity model. The sector- 
selective transitions studied here are then simply different 
instantiations of the general phenomenon of the "orbital- 
selective Mott transition," extensively discussed in the 
context of transition-metal oxide a^^'^°i^-^'^^i'^^'^°i^-^ with 
partially filled d levels. This connection was explored by 



Ferrero and co-workersi^iH and also by Liebsch^i and 
provides both insights into the physics of the cluster dy- 
namical mean-field theory and insights back to the single- 
site physics of multiorbital transition-metal oxides. 



However, we note that when expressed in the orbital 
basis the interactions in the cluster studied here are con- 
siderably more complicated than those normally studied 
in the orbitally degenerate transition-metal oxide case. 
Liebsch et al. found that the four-site clusters they stud- 
ied did not exhibit momentum-selective transitions. The 
presence of two interaction-driven transitions at half fill- 
ing seems (as noted above) to follow naturally from the 
different bandwidths of the two relevant sectors while the 
persistence of the orbital-selective effect on doping is sim- 
ilar to the two-orbital results presented in Ref. [50- We 
note, however, that neither the second-order nature of 
the first transition nor the collapse of the two transitions 
into one as the t' parameter is increased seem to have 
been noticed in the orbital-selective Mott-transition lit- 
erature. These issues seem worth exploring in the latter 
context. Further, the essential role apparently played in 
our calculations by the pinning of the density to 1/2 in 
the insulating orbitals requires further investigation. In 
the two-site calculations of Ferrero et al. this pinning 
was not found. 



From the point of view of the orbital-selective Mott 
transition the models studied here have a special feature, 
which is that the different "orbital" (A') sectors are de- 
coupled both by the symmetry of the impurity model 
and by the self-consistency condition. In other cluster 
DMFT implementations (for example, cellular dynamical 
mean-field theory) this may not be the case. In the con- 
text of the orbital-selective Mott phenomenon this decou- 
pling has been shown to have an important consequence; 
it rules out a Kondo coupling between the gapped and 
the ungapped orbital which, if antiferromagnetic, might 
destabilize the orbital-selective phase at temperatures be- 
low the Kondo scaleJ^ It is therefore natural to ask if the 
orbital-selective phase we find survives beyond the DCA 
approximation we have used here. We note that the tran- 
sitions we find in the cluster DMFT context have an ad- 
ditional feature not yet explored in the orbital-selective 
literature the orbital which first becomes insulating (sec- 
tor C) is two fold degenerate. The insulating orbitals 
may therefore have a more complex spin structure than 
in the non-degenerate case, and, in particular, it seems 
likely that the singlet state of the two orbitals will be fa- 
vored, in which case a Kondo effect would not necessarily 
occur. Cellular dynamical mean-field theory calculations 
(in which the self-consistency condition couples the sec- 
tors) would be very valuable, and it may be interesting 
to explore the sector degeneracy effects further in the 
selective Mott-transition literature. 
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IX. SUMMARY 

Building on previous rcsultsi^ we have shown in 
this paper that the eight-site DCA approximation to 
the intermediate-couphng Hubbard model yields results 
which exhibit a remarkable qualitative correspondence 
to the pseudogap effects occurring in high-temperature 
superconductors. For parameters relevant to high- Tc su- 
perconductors we find that hole doping but not electron 
doping leads to a two-stage metal-insulator transition in 
which a gap opens first around the (tt, 0) point. Our con- 
clusions are based on an analysis of directly measured 
quantities and do not involve analytical continuation or 
interpolation of self-energies. In this concluding section 
we summarize our results and place them in the context 
of other studies of high-temperature superconductivity 
and of other dynamical mean-field analyses. This section 
includes pointers back to relevant sections of the text and 
may be read independently of the rest of the paper. 

A crucial finding of Ref.lTslwas that for couphng strong 
enough that the model is insulating at n = 1 a two- 
stage metal-insulator transition occurs as the insulating 
phase is approached from the regime of large hole dop- 
ing. As the doping is reduced below a critical doping 
Xc (~ 0.10 for U = 7t and increasing weakly as U is 
increased) a gap appears in the momentum sector con- 
taining the "antinodal" region (0, tt) while the states in 
the nodal region remain gapless until the half-filled point 
is reached. The resolution available to us does not per- 
mit us to distinguish between a full gap (density of states 
= for a nonzero range of energies) and a soft or pseudo- 
gap where the density of states vanishes only at one fre- 
quency or is strongly suppressed over a range of frequen- 
cies but does not vanish. We found that in this "sector- 
selective" region the antinodal gap is small in comparison 
to the correlation gap found at half filling, although it 
does grow rapidly as doping is decreased. The small gap 
magnitude but rapid increase can be seen from the rel- 
atively weak temperature dependence of the many-body 
density of states (Fig. [3]) and from the Matsubara-axis 
self-energies fFigs. [T^ and [T^ . For an electronic disper- 
sion with particle-hole symmetry the two-stage transition 
occurs for both electron and hole doping but for particle- 
hole asymmetries of the magnitude believed^ to be rele- 
vant for high-temperature superconductors we find that 
the two-stage transition occurs only on the hole-doping 
side (Fig. [7]). These two findings are the essence of the 
high- Tc pseudogap phenomenon and thus our results add 
support to the widely held belief that the intermediate- 
coupling Hubbard model contains much if not all of the 
essential physics of high- Tc superconductivity. 

The "antinodal gap" phase appears to be non-Fermi- 
liquid in two senses; of course the gapping of the antin- 
odal sector implies that the Fermi surface violates the 
Luttinger theorem associated with the full carrier den- 
sity. Further, ungapped states in the nodal sector are 
non- Fermi-liquid in the sense that the many-body density 
of states appears (Fig. [3]) to assume a low-temperature 



value which is neither the non-interacting value (as ex- 
pected if the states in this sector were Fermi-liquid like) 
nor vanishing (as expected if the sector were gapped) and 
in the sense that the imaginary part of the Matsubara- 
axis self-energy appears not to vanish as iujn — *■ at the 
temperatures accessible to us fFig. [T3)) . 

The antinodal gap regime is separated from a large- 
doping Fermi-liquid phase by a zero-temperature transi- 
tion whose existence is inferred from a crossing of the 
temperature dependence of the many-body density of 
states. The transition is apparently continuous, char- 
acterized by a gap which opens smoothly as the doping 
is decreased below the apparent quantum critical point 
fFigs. [9] and [T2)) . Because the transition does not involve 
an order parameter it is a smooth crossover at T > 0. 
An apparently rather similar transition occurs if the in- 
teraction is varied at density n = I and the particle-hole 
symmetry breaking is not too large (Figs. [8] and [9]). 

The second transition, at which the nodal sector devel- 
ops a gap, is first order (characterized by discontinuities 
in the gap amplitude) in the interaction-driven case. In 
the doping-driven case the behavior of physical quantities 
appears smooth but the chemical potential is tuned into a 
large-amplitude pre-existing gap in the nodal sector (see 
Fig.IHD. 

An antinodal gap which varies smoothly with doping 
and vanishes above a critical doping level, is an essential 
feature of the pseudogap observed in hole-doped high- 
Tc materials^. Interestingly, the many-body density of 
states obtained by analytical continuation of quantum 
Monte Carlo data for four-site clusters (which do not 
provide the momentum resolution needed to distinguish 
the nodal and antinodal sectors) revealed in the doped 
case a small gap pinned to the Fermi levelJ^ Ferrero et 
ali^'^ subsequently found a similar feature in a two-site 
calculation. We speculate that this is the expression, in 
the smaller clusters, of the two-stage transition which 
is clearly revealed in the larger eight-site cluster. Fur- 
ther investigation of the relation between two-, four-, and 
eight-site cluster results is in progress. 

The transition at which sector C becomes gapped ap- 
pears to be different from the quantum phase transition 
identified by Vidhyadhiraja et alr^ and the related phe- 
nomena studied by Liebscbi^ using a four-site cluster. 
The Vidhyadhiraja transition occurs at the higher doping 
X ~ 0.15 (corresponding to chemical potential /i w ~2t) 
where within our resolution no effects occur. However, 
we note that the transition identified by Vidhyadhiraja 
et al. was obtained for a different cluster geometry and 
more importantly involves a change in the self-energy 
from a Fermi-liquid to a marginal Fermi-liquid form; such 
a change would lead only to subtle effects in the quan- 
tities we are able to measure accurately. The range of 
temperatures we can access may also be too small. 

It is important to note that the transition is not driven 
by the van-Hove singularity in the density of states. As 
can be seen from Fig. [TOl at t' = — 0.15i the transition 
to the antinodal gap phase occurs after the rcnormalizcd 
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Fermi level has passed through the Van Hove point. A 
generic feature of weak-coupling theories based on per- 
turbation about a Fermi-liquid phase is the importance 
of Van Hove physics. The absence of a clear association 
with the Van Hove point identifies the transition we find 
as an intermediate or strong-coupling phenomenon. 

The sector-selective transitions that we find are as- 
sociated with a pinning of the sector occupancy to the 
half- filled value risoctor = 1/2. The pinning is a T ^ 
effect (consistent with the lack of an order parameter 
and nonexistence of a T > transition) and because of 
the small gaps the approach to the T ~ limit can be 
slow. The pinning of the density and the slow approach 
to the r = value can be seen in Fig. 21 The importance 
of the pinning of the density is also revealed by the de- 
pendence of our findings on the particle-hole symmetry 
breaking parameter t' . The sector-selective phenomenon 
occurs both as the interaction strength is varied at den- 
sity n = I (Fig. and as the density is varied at strong 
interaction (Fig. [7]). 

In the interaction-driven case a two-stage transition 
occurs only for sufficiently small. At t' ~ and 
n = 1 both the nodal and antinodal sectors are half filled 
at all values of the interaction strength and the smaller 
bandwidth of the antinodal sector (Fig. [5]) explains why 
this sector becomes insulating at a smaller interaction 
strength than does the nodal sector. However, for t' 
in the small U fully metallic phase both sectors have oc- 
cupancy different from the half-filled value while in the 
sector-selective phase the gapped sector (and therefore 
also the ungapped one) have density equal to the half- 
filled value. As U is increased , momentum differentia- 
tion (a difference in electron self-energy between sectors) 
causes charge to be transferred between sectors. Thus 
the sector-selective phase is energetically disfavored by 
the need to overcome the kinetic energy and equalize the 
charge distribution between sectors. At small \t'\ this 
effect leads only to a weak increase in the critical inter- 
actions strength but as \t'\ is increased the critical U re- 
quired to equalize the charge distribution becomes larger 
than the value needed to open a gap in the nodal sector, 
leading to one first-order transition. 

A very similar phenomenon may be seen in the doping- 
driven transition (Fig. [7]). For small \t'\ sector-selective 
transitions occur for both electron and hole-doping but 
as t' /t becomes more negative (positive), the transitions 
on the electron-doped (hole-doped) side merge and be- 
come first order, again because the density difference be- 
tween sectors becomes too large to sustain without gap- 
ping both sectors. This physics is revealed most clearly in 
Fig. m where for t' /t = —0.3 we see that as the chemical 
potential is raised across the insulating gap the density 
of the B sector actually drops below the half-filling value 
even though the mean electron density increases. Simi- 
larly, on the hole-doped side the sector-selective transi- 
tion occurs when the chemical potential reaches the value 
at which the antinodal sector becomes half filled. 

The close association between commensurate filling 



and the antinodal gap leads us to identify the antinodal 
transition as being of the Mott type. We also observe in 
passing that no evidence for "nematic" density fluctua- 
tions or a Pomeranchuk-type Fermi surface instability is 
found (Fig. llSp . If the transition were driven exclusively 
by short-ranged antifcrromagnctic order there would, as 
far as we are aware, be no reason for the sector density 
to be pinned. We observe, however, that the onset of 
the sector-selective phase is associated with a strong in- 
crease in quasi thermal, long correlation length antiferro- 
magnctic spin fluctuations (Fig. [H)) . As in the single-site 
dynamical mean-field theory, the onset of a Mott phase 
is associated with local-moment formation and the pres- 
ence of the local moments will of course lead to strongly 
enhanced spin correlations. Interestingly, we see no ob- 
vious sign of a spin gap. While one expects the eight-site 
cluster to have a low-spin ground state, the energy dif- 
ference between lowest-lying singlet and triplet states is 
evidently below our temperature resolution. 

In summary, we have presented evidence that 
the intermediate-coupling Hubbard model exhibits a 
momentum-selective, multistage Mott transition leading 
to a pseudogap phenomenon in remarkable qualitative 
agreement with high-Tc materials. The systematics of 
our findings strongly suggest that a weak-coupling in- 
terpretation is not tenable. Important open questions, 
presently under investigation, include the relation of the 
present eight-site cluster results to those obtained on 
smaller and on larger clusters. In this context we note 
that studies of non particle-hole symmetric clusters (ei- 
ther doped or with t' ^ 0) encounter a fermion sign prob- 
lem which becomes progressively more severe for lower 
temperatures, stronger interactions, and larger clusters 
so unless some remarkable new methodological innova- 
tion appears extension of our results too much beyond 
the parameters studied here will not be possible. How- 
ever, working with presently accessible cluster sizes and 
interaction strengths, an extension of our results to in- 
clude long-ranged magnetic and superconducting order 
is possible and is an important open question. 

Recently, we became aware of work of Khatami et al.^ 
which reported DCA calculations of the compressibility 
dn/dfi indicating a sharp first-order transition with asso- 
ciated hysteresis at a small electron doping for |<'| > 0.2. 
We associate this transition with our findings (Fig. 2] and 
[7]) that for \t'\ > 0.2 and large U the doping-driven metal- 
insulator transition is first order on the electron doped 
side. We expect that this first-order transition will lead 
to hysteresis and anomalous behavior of dn/dfi at low 
dopings, on the order of the jump in density shown in 

Fig.m ' 
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